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Abstract. We are concerned with multidimensional stochastic balance laws. 
We identify a class of nonlinear balance laws for which uniform spatial BV 
bounds for vanishing viscosity approximations can be achieved. Moreover, we 
establish temporal equicontinuity in of the approximations, uniformly in 
the viscosity coefficient. Using these estimates, we supply a multidimensional 
existence theory of stochastic entropy solutions. In addition, we establish 
an error estimate for the stochastic viscosity method, as well as an explicit 
estimate for the continuous dependence of stochastic entropy solutions on the 
flux and random source functions. Various further generalizations of the results 
are discussed. 



1. Introduction 

We are concerned with the well-posedness and continuous dependence estimates 
for the stochastic balance laws 

dtu{t,K) + \/ ■{{u{t,xj)^a{u{t,x))dtW{t), xeR'^,t>0, (1.1) 

with initial data: 

u(0,x) = ■uo(x), xeR''. (1.2) 
We denote by V and A the spatial gradient and Laplacian, respectively. 

Equation (jl.ip is a conservation law perturbed by a random force driven by a 
Brownian motion W{t) — W{t,uj), uj £ V,, over a stochastic basis (fi, J^, {J^t}t>o ' P)^ 
where P is a probability measure, is a cr-algebra, and {J^t}(>Q is a right-continuous 
filtration on {^,J^) such that contains all the P-negligible subsets. 

The initial function uo(x) is assumed to be a random variable satisfying 

^ [|I"oIIlp(r<') + l^olsy(R<i)] < p=l,2,---. (1.3) 

Regarding the flux f = (/i, • • • , /d) : M ^ M"*, we assume that /; e C^{R), 
i ~ 1, . . . , d, and that each fi has at most polynomial growth in u, i.e., 

|/*(")l <C{1 + l^r) for some finite integer r > 0. (1.4) 

In this paper we focus mainly on the class of noise functions a for which there 
exists a constant C > such that 

cr(0) = 0, \(7{u) - (7{v)\ <C\u-v\ yu,v£R. (1.5) 

This can be generalized to wider classes for different results in terms of existence, 
stability, and continuous dependence, respectively; see Section [5] for more details. 
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One reason for requiring cr(0) = is that it follows from the L^-contraction principle 
that E[\\u{t, is finite. Similarly, the Lipschitz continuity of cr(u) is required 

for the existence and uniform estimates of solutions. 

Stochastic partial differential equations arise in a number of problems concerning 
random-phenomena occurring in biology, physics, engineering, and economics. In 
recent years, there has been an increased interest in studying the effect of stochastic 
forcing on solutions of nonlinear stochastic partial differential equations. Of specific 
interest is the effect of noise on discontinuous waves, since these are often the 
relevant solutions; an issue of particular importance concerns the well-posedness 
(existence, uniqueness, and stability) of discontinuous solutions. 

The fundamental fluid dynamics models are based on the compressible Navier- 
Stokes equations and Euler equations. However, abundant experimental observa- 
tions suggest that the chaotic nature of many high- velocity fluid dynamics phenom- 
ena calls for their stochastic formulation. Indeed, in these flows with large Reynolds 
numbers, microscopic perturbations get amplified to macroscopic scales giving rise 
to unsteady fiow patterns that deviate significantly from those predicted by the 
classical Navier-Stokes/Euler models, and more viable models seem to be the sto- 
chastic Euler or Navier-Stokes equations. In the present paper we are interested in 
nonlinear hyperbolic equations with stochastic forcing, so-called stochastic balance 
laws. These balance laws can be viewed as a simple caricature of the stochastic 
Euler equations. 

Some efforts have been made in the analysis of nonlinear stochastic balance laws. 
When cr = 0, becomes a nonlinear conservation law for which the maximum 

principle holds. A satisfactory well-posedness theory is now available (cf. [5j). 
In [TU], a one-dimensional stochastic balance law was analyzed for uq in L°° and 
compactly supported a — a{u), which ensures an L°° bound. A splitting method 
was used to construct approximate solutions, and it was shown that a subsequence 
of these approximations converges to a (possible non-unique) weak solution. 

For general a, the maximum principle is no longer valid. Indeed, even for L°° 
initial data uq, the solution is no longer in L°° generically. For cr = (j(t,x) in 
Ct(W^'°°) and with compact support in x, Kim [12] established the existence and 
uniqueness of entropy solutions in the one-dimensional case; see also |22j . For more 
general a — (j{x, u) depending on u and for multidimensional equations in the U' 
framework, the uniqueness of strong stochastic entropy solutions was first estab- 
lished in Feng-Nualart [5] , but the existence result was restricted to one dimension; 
see the recent paper Debussche-Vovelle '6' for multidimensional results via a kinetic 
formulatiorQ. For the theory of deterministic conservations laws, see |21) . 

One of our main observations is that uniform spatial BV bounds are preserved 
for stochastic balance laws with noise functions <j{u) satisfying (jl.Sp . This yields 
the existence of strong stochastic entropy solutions in n BV , as well as in L^, 
for multidimensional balance laws (11.11) . Furthermore, we develop a "continuous 
dependence" theory for stochastic entropy solutions in BV ^ which can be used, for 
example, to derive an error estimate for the vanishing viscosity method. Whenever 
(T = c(x, u) has a dependency on the spatial position x, BV estimates are no longer 
available, but we show that the continuous dependence framework can be used to 
derive local fractional BV estimates, which in turn can be used, as before via a 
temporal equicontinuity estimate, to establish a multidimensional existence result. 



We became aware of this paper after our main results were obtained. 
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Besides providing an existence result in a multidimensional context by standard 
methods, one reason for singling out the class of nonlinear balance laws defined by 
(|1.5p is that it makes a natural test bed for numerical analysis, without having to 
account for all the added technical complications in a pure framework. Moreover, 
by assuming cr(a) = a(b) = for some constants a < one ensures that the solution 
remains bounded between a and b if the initial function uq does so. Consequently, it 
is possible to identify a class of stochastic balance laws for which n BV, or even 
L°° n BV, supplies a relevant and technically simple functional setting, tailored for 
the construction and analysis of numerical methods. 

For other related results, we refer to Sinai [19 and E-Khanin-Mazel-Sinai [7] 
for the existence, uniqueness, and weak convergence of invariant measures for the 
one-dimensional Burgers equation with stochastic forcing which is periodic in x, 
as well as the structure and regularity properties of the solutions that live on the 
support of this measure. We also refer to Lions-Souganidis (TB] for Hamilton- Jacobi 
equations with stochastic forcing and the so-called "stochastic" viscosity solutions. 

We employ the vanishing viscosity method to establish the existence of stochastic 
entropy solutions. To this end, consider the stochastic viscous conservation law 

dtu%t, x) + V • f x)) = a{u%t, x.))dtW(t) + eAu^t, x) (1.6) 

for any fixed £ > 0, with initial data 

w^(0,x) =u^(x), xeR'', (1.7) 

where 'Uq(x) is a standard mollifying smooth approximation to mo(x) with 



E 



< E 



|wo(x)|Pc;x 



and, if uo e BV(M'^), 



E [ |Vu^(x)| dx <E j |Vuo(x)| dx . 

In addition, E [J^^ |V^Uq(x)| dx] < oo, i.e., |V^Uq| is integrable for each fixed e. 

With regard to the viscous equation (|1.6p . we should replace {f,(j) by appro- 
priate smooth approximations (f^,(T'^). However, mainly to ease the presentation 
throughout this paper, we will not do that but instead simply assume that (f , a) 
are sufficiently smooth (cf. [9]) in order to ensure the validity of our calculations. 
At times, we will do the same with the initial data. 

The existence of global smooth solutions to (|1.6|) - (|1.7|) is established in along 
with the following uniform estimates for p > 1 and T > 0: 



sup sup 

e>0 0<t<T 



£^[ll"'(*r)lli.(K.) 



1 +sup£; 


'^1 


^ e>0 


Jo 



|Vu^(t,-)|li.( 



< oo. 



(1.8) 



The solution satisfies 



u%t,Ti) = I Ge{t,x-y)ua{y)dy 

f f G,{t-s,x-y)V ■i{u'{t,y))dyds 
Jo Jw 

Ge(i - 5, X - y)<j{u'{s, y)) dy dW{s), 



(1.9) 
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where Ge(<,x) is the heat kcrnch 

Using ([Q]) and ([L8l) -(fr9 l) . it follows that, for each fixed £ > 0, 

•)] < CXI for any finite T > 0, 



E[\\iV,AX 



(1.10) 



\L^{(0,T)xR'^)\ 

that is, Vu*^ and V^w^ are integrable for each fixed e > 0. 

With different methods, we will later prove an e-uniform spatial BV estimate. 

The remaining part of this paper is organized as follows: In Section [2l we prove 
the uniform spatial BV bound for stochastic viscous solutions u^{t,x). Based on 
the BV bound, we establish the equicontinuity of x) in i > 0, uniformly in the 
viscosity coefficient £ > 0, in Section[3l With these uniform estimates, we establish 
the existence of stochastic entropy solutions in n BV, as the vanishing viscosity 
limits for problem (ll.6l) - (|1.7l) with initial data in L^OBV, in Section|4l Combining 
this existence result with the i^-stability theory in Feng-Nualart [9] leads to the 
well-posedness in for problem (|l.ip - (|1.2p . We further establish estimates for the 
"continuous dependence on the nonlinearities" for BV stochastic entropy solutions 
in Section [5l which also leads to an error estimate for (ll.6p - (|1.7p . Various further 
generalizations of the results are discussed in Section |6l 

2. Uniform spatial _BT^-estimates 

As indicated in Section [1] we have known the regularity and the uniform L^- 
estimate (|1.8I) {p > 1) for the viscous solutions u'^{t, x) of (|1.6p - (ll.7p . In this section, 
we establish the uniform L^-estimate for Vu^, that is, the uniform BU-estimate of 
u^{t,x.) in the spatial variables x. 

Before we do that, let us indicate why BV estimates do not seem to be available 
when the noise coefficient function a = a{x,u) depends on the spatial position x, 
even if that dependence is C°° (see Section [5] for fractional BV estimates). To this 
end, it suffices to consider the simple stochastic differential equation: 



du — a{x,u) dW{t), u{0) — uq{x), 



X e 



where we have dropped nonlinear transport effects and restricted to one spatial 
dimension. The spatial derivative v — dxU satisfies 

dv = {<Ju(x, u)v + ax{x, u)) dWli:,). 

Let 77 be a C^-function. By Ito's formula, 

1 2 
dr]{v) = r]'{v)(au{x, u)v + a^ix, u)) dW{t) + --q" (v)(au(x, u)v + ax{x,u)) dt. 

Integrating in x and taking expectations, it follows that 



E 



?7(u(i)) dx 



E 



E 



r/(u(0)) dx 
1 



rj" {v)(^au{x,u)v + ax{x,u)) dxds 



Modulo an approximation argument, we can take ry(-) as |-|. Unless ax = 0, the 
second term on the right-hand side does not seem to be controllable (this term 
vanishes when ax =0). 



STOCHASTIC BALANCE LAWS 



5 



Let US now continue with the derivation of the BV estimate for ()1.6p . We will 
need a -approximation of the Kruzkov entropy. Let 77 : M — > R be a C^-function 
satisfying 

7?(0) = 0, ?7(-r) = J7(r), fj' hr) = -fj' (r), 7?" > 0, 

and 

-1, when r < —1, 

e [-1, 1], when \r\ < 1, 
-1, when r > 1. 

For any p > 0, define the function r/s : M — > M by 



fj'ir) 



(2.1) 
(2.2) 



Vpir) = Pvi-)- 
P 



Then 



where 



M2 



|r| - Mip < rjp{r) < \r\ , \Tj';{r)\ < ^1|,.|< 
Ml = sup I \r\ - f]{r)\, = sup \fj"{r)\ 



|r|<l 



r|<l 



(2.3) 

(2.4) 
(2.5) 



We will frequently utilize the Burkholder-Davis-Gundy inequality, which we now 
recall. For p > 0, there exists a constant C = Cp such that, if Mt is a continuous 
martingale and t a stopping time, then 



E 



sup|A/,|f 



< 



CpE [(Af)^ 



p/2 



where {M)t is the quadratic variation of Mt. 

Theorem 2.1 (Spatial BV estimate). Suppose that (|1.3p - (|1.5p hold. Let u^(i,x) 
be the solution of (jl.6p - (|1.7p . Then, for t > 0, 



E 



/ |Vu^(t,x)|dx 



< E 



|Vmo(x)| dx 



< l; 



|Vuo(x)| dx 



Proof. Taking the derivative of (|1.6p with respect to Xi, 1 < « < d, we obtain 
a*«)+V. (f'(7.^(t,x))0 =a'(7.-(t,x)XaiM^(<)+£A«J. 
Applying Ito's formula to rip{u%.) yields 



+ 77;«J(£A<-V.(f'(7.>^J) 

We observe that 

£r7;«)A«J = £ (V • (i«)V<J - <«J |V<|') 
= e(Ar/p«)-<«)|V<j2) 

< £A7;p«J, 



(2.6) 



(2.7) 



by using the convexity of r^p, and interpreting A77p(m| .) in the distributional sense. 
Here we have used that < ? < d, are integrable (cf. (jl.lOp ) so that they 

vanish at infinity, which leads to the vanishing boundary terms in (j2.7p . 
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Integrating ()2.6p with respect to x, using ()1.10|) and (|2.7p . and noting that 



r;'«Jcr'(u^X^ dWis)dx 



IR'' Jo 

is a martingale, we arrive at 



E 



E 



?7p«^(0,x))fix 



< E 



1 



i«)V-(f'KKJdxds 



(2.8) 



Now we send p — > in (|2.8p . By the dominated convergence theorem, 



< E 



/ |u|.(t,x)| dx 
/ K.(0,x)| dx 



Um E 



77;«)V-(f'(^.^)Odxds 



hm — i? 

P^o 2 



JR'' 



/ |<,(0,x)| dx 



+ /l +/2. 



For the Ii term, 



= hm 

p— !-0 



E 



hm 

p^O 



< C hm 

p^O 



n'7p«.X.V<_ •f'(w'')dxds 



"'Rt' 



Notice that 



and 



\^xi \ -X[-p,p]{Uxi) ^ for a.e. (i,x) almost surely as p — ?► 0, 



-X[-p,p]«)|V<||f'(u^)| 



< C Vu 



where the right-side term of the inequality is integrable and independent of p > 0. 
Then the dominated convergence theorem implies that =0. 

Next we consider 12- By condition (|1.5I) and estimate (|2.4I) . we have 
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On the other hand, since | is integrable and independent of p > and 

I I ■'■||tt= |<p} ~^ ^ (^'■''■) almost surely as p — > 0, 

the dominated convergence theorem again implies I/2I — 0. This concludes the 
proof. □ 

3. Uniform Temporal L^-Continuity 

In this section, we establish the uniform temporal L^-continuity of w'^(t,x), 
independent of the viscosity coefficient e > 0. 

Theorem 3.1 (Temporal L^-continuity). Suppose that (|1.3l) - (jl.5p hold. Letu'^{t,yL) 
be the solution of (jl.6p - (11.71) . Let D C M'' be a bounded domain in M'' and T > 
finite. Then, for any small At > 0, there exists a constant C > independent of 
At .such that 



E 



" i-T-At I- 

/ / |u''(i + At,x) - w''(t,x)| dxdt 

"'0 JD 



T~At 
"T) 

< C{Aty/^ -> as At -> 0. (3.1) 



Proof FixAt>0. Fori e [0,T- At], set w^(t,-) := u"" (t + At, ■) - {t, ■). Then, 
for any ip G L°°(0, T; Cg°{D)), we have 



/ x)(y5(i, x) dx 

JD 

= j (^j dsu'^{s,x)ds^(p{t,x) d:!C 

pt+At P 

= / / f(u''(s,x)) • V<^(t,x)dxds 

rt+At p 

— e / / Vw'^(s, x) • Viy9(t, x) dxds 

Jt JD 



(3.2) 



t+At 

/ a{u^s,:ic))(p{t,x.)dyidW{s). 
t JD 

For each te [0,T - At], take S > 0, set D^s {x e L» : dist(x, dD) > S}, and 
denote by XD^si') its characteristic function. 

Let J e C;?°(R'') be the standard mollifier defined by 

J(x).f^^P(^) ^^l-Kl' (3.3) 
lO if |x| > 1, 

where the constant C > is chosen so that J^^^ J{x.)dx ~ 1. For each S > 0, we 
take 

-d / 7/x:iy 



if -.^ ipsit,x.) = S-"" / J(^i^)sgn(u;(t,y))xD_.(y)ciy 

JR'' 

in p.2p . It is clear that ||</'<5|lL°°(r>) + ^ II^'/'5|Il°°(_d) — C*, uniformly in t, for some 
constant C > independent of 5 > 0. 
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Integrating p.2p in t from to T — At yields 



T-At p 

/ |w^(t,x)| dxdt 

T-At ^t+At 



/ / / f(u^(s,x)) • V<y95(i,x)dxdsdt 

^T-At pt+At p 

/ / / eVu^(s,x) • Vi^^(<,x) dxdsdt 

Jo Jt Jd 

f7(u^(s,x))(p5(i,x)dx)dT/F(s) ) 
w^{t, x)(w'^(i, x) — (^^(t, x)) dxdt 



Jt J D 

T~At / ft+At 



\Jt ^ JD 

T-At 



JD 



We examine these parts separately. 

Thanks to the polynomial growth of f and (|1.8p . 



E [H] 



<G^m\L^(D.i,,T))<C{T,D)^. 



For the term we have 









E [4] 




irw. 







X E 



T-At 



< CAt E 





T-At 







pt+At 

y/e\\/u''{s,x)\ds ] dxdt 
e|V(p^p dxds 
\\/(ps\'^dx ds 



<C{T,D)^, 



where the second inequality follows from the energy estimate (11.8 



supi? 

£>0 



: / |lVw^(t,x)||^2(R.) 
Jo 



dt 



< oo. 
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For the term /|, by the Burkholder-Davis-Gundy inequality apphed to the mar- 
tingale < At h^- Jj*^^* ( Jjj a{u^ {s,x))(ps{t,-x.) dx) dW{s), we have 



E [4] 



pT-At 

<C E 
Jo 



t+At ^ j- N 2 \ 2 

C7(u'^(s, x))i^5(t, x) dx ) ds\ 



<C\E 



t ^ JD 

-T-Al pt+At 



dt 



<C[ E 



Jt J D 

At rT-At 



< Cy/Ail E 



< C\fKt ( E 

< cVm. 



a{u^{s,x)(ps{t,yi)) dxdsdt 
II I ((7(u^(s + i,x)) dxdtds 

Jo Jo J D 

f f {a{u%t,x)fdxdt 
Jo Jd 

C [ \u%t,x)\ 
Jd 



'dxdt 



where we have used that sup^>o-^ [ll^^lOlli] < uniformly in t > 0. 
This L^-bound also implies 



E 



fl 

Jo Jd\d 



\u%t,x)\dxdt 



<C {E[\\u^\\l]y [E[ / dxdt] 
\ Jo Jd\D_2s 

< C\/6. 
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Hence, 
E[ll] 

< 2E 



+ E 


[L 








< CV6 








+ E 


'j 








n 


< cVs 


+ 


CE 



T-At 



l-T-At 



\w{t, x)| dx.dt 



\w{t,^)\-w{t,^) / rt(^)sgn(ti;(i,y)) 



dy dx dt 



D^2 



T-At 



S-'^Ji^) |w(t,x)| -u;(i,x)sgn(w(i,y))|dydxdi 



D-2S 



5-''j{^)\w{t,^) - w{t,y)\dy d^dt 



< cVd + CE 



I J{z) [ f \u%t,x) - u%t,x- Sz)\dxdtdz 



< cs^ +4S< cVs, 

where the third inequahty fohows from ||a| — asgn(&)| < 2|a— b\ for any a,b 
The fifth inequahty fohows, since u'^ belongs to BV in x. 

Setting p{At) = Ms>o {Ci^ + C2(At)^ + €36^, it fohows that 



T-At 



|w(t,x)| dxdt < p(At). 



D 



The function p(-) reaches the infimuni at S = C{At) 3 ^ and hence 

fT-At 

\w{t,x)\dxdt < C{At)i ^ as At ^- 0. 



D 



□ 



Remark 3.1. Since Brownian sample paths are a- Holder continuous for every a < i, 
a fractional order in the temporal -continuity in p.ip is expected. The proof of 
Theorem 13. II uses an idea due to Kruzkov [131 . 



4. Well-Posedness Theory in Lp 

Before we introduce the relevant notions of generalized solutions, let us define 
what is meant by an entropy-entropy flux pair (77, q), or more simply an entropy 
pair, namely a function 7y : R — >■ M such that 77',??" have at most polynomial 
growth, with corresponding entropy flux q defined by q'(w) — r]'{u)t'{u). An 
entropy pair is called convex if if' {u) > 0. 

Definition 4.1 (Stochastic entropy solutions). A {J"(}j^p-adapted, L^(R'^)-valued 
stochastic process u = u(t, x; to) is a stochastic entropy solution of the balance law 
p.ip with initial data (|1.2p provided that the following conditions hold: 
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(i) forp=l,2,..., 



0<t<T 



sup E[\\u{t)\\l^^^a^] <oo, foranyT>0; 



(ii) for any convex entropy pair {t], q) and any < s <t, 





for all if e C~(]R'^), > 0, where (• • • ) dW{T) is an Ito integral. 

To motivate the next definition, let us make a formal attempt to derive the L^- 
contraction property for stochastic entropy solutions. To this end, consider smooth 
(in x) solutions to the onc-dimcnsional problems: 



Let ?7(-) be an entropy. An application of the chain rule (Ito's formula) now yields 



where the last term is a martingale. Choosing ry(-) = |-| yields r]"{-) — do and the 
two "jj" terms" vanish. Consequently, after integrating and taking expectations, 
we arrive at the X^-contraction (conservation) principle: 



Of course, for non-smooth solutions, the Ito formula is not available and we 
should instead derive the -contraction principle from the (stochastic) entropy 
inequalities via Kruzkov's method. 

Attempting precisely that, we write the entropy condition for u{t) = u{t,x;uj) 
with the entropy rj{u{t) — v{s,y;Lo)), where ?;(s,?y;cj) is being treated as a constant 
with respect to {t,x). Similarly, write the entropy condition for v{s) = v{s,y;oj) 
for the entropy r]{v{s) — u{t,x;co)), with u{t,x;uj) being constant with respect to 



du + dxf{u) dt = a(u) dW, w|t=o = uq, 

dv + d.,f{v) dt = a{v) dW, v\t=o = Vo- 



Subtracting the two stochastic conservation laws yields 

d{u -v) = - [d4f{u) - f{v))] dt + [a{u) - a{v)] dW. 



dr]{u-v)= -d,{rj'{u-v){f{u)-f{v))) 



+ rj"{u-v){f{u)-f{v))d,{u~v) 

1 2 

+ 2 V(w - v) {a{u) - a{v)) dt 
+ ri'{u - v) (<t(m) - a{v)) dW, 
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(s,y). Take ri{-) — |-|, and then q{u,v) — sgn(u — v){f{u) — f{v)). After adding 
together the two entropy inequahties, we formally obtain 

{dt + ds)\u - v\ 



< 



dt ds 



{d,+dy){sgn{u-v){f{u)- f{v))) 

+ sgn(u(i, x) — v{s, y))(j{u{t, x)) dW[t) ds 
- sgn(u(i, x) - v{s, y))a{v{s, y)) dW{s) dt. 



Depending on i < s or i > s, one of the last two terms are not adapted, and this 
causes a problem for the Ito integral. In particular, by taking the expectation of 
the above inequality, only one of the last two terms vanishes. Moreover, to write 

in the favorable form: 

1 2 

-S{u - v) {(t{u) - (t{v)) , 

we are missing the cross term 2a{u)a{v). These difficulties can be effectively han- 
dled by the notion of "strong" stochastic entropy solutions. 

Definition 4.2 (Strong stochastic entropy solutions). An {J^t}j^Q-adapted, L^(M'^)- 
valued stochastic process u — u{t) — u(t, x; uj) is a strong stochastic entropy solution 
of the balance law (|l.ip with initial data (|1.2p provided m is a stochastic entropy 
solution, and the following additional condition holds: 

(iii) for each { J'tlj^Q-adapted, L^(R)-valued stochastic process u — u{t) — 
'u(i,x;a;) satisfying 



sup E \\\u{t)\\ 

0<t<T 



< oo for any T > 0, p = 1, 2, 



and for each entropy function S : 



with 



S{r;v,y) 



S'{u{r, x) — v)a{u{r, x))(^(x, y) dx. 



where r > 0, v € R, y e R'^, and ip e C^{R'^ x R'^), there exists a 
deterministic function A{s,t), < s <t, such that 



E 



< E 



S{t;v = u{t, y),y) dW{T)dy 



dyS{T; V = u{t, y), y)a{u{T, y)) dy dr 



A(s,t), 



where A(-, •) is such that, for each T > 0, there exists a partition {ti}^^ 
of [0, T], = to < ti < • • • < t„, = T, so that 



lim VA(t„t,+i) 

ma.x.\ti^i—ti \ 



0. 



The notion of strong stochastic entropy solutions is due to Feng-Nualart [5] , who 
proved the L^-contraction property for these solutions: 

E \\u(t) - f < E [||«o - t'olUiCR.^)] for t > 0, (4.1) 
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where u{t) is any stochastic entropy solution with u\t=o — uo and v{t) is any strong 
stochastic entropy solution with v\t=o — vq. In (|4.ip . the entropy |-| can be replaced 
by (■)^, yielding the L^-comparison principle. 

Feng-Nualart |9] employed the compensated compactness method to prove an 
existence result in the one-dimensional context. The following theorem provides 
the existence of strong stochastic entropy solutions for a class of multidimensional 
equations. 

Theorem 4.1 (Existence in LP D BV). Suppose that (fOI) -(fTS ]) hold. Then there 
exists a strong stochastic entropy solution u of the balance law (jl.ip with initial 
data (jl.2p . satisfying 

-E' •)li3y(R'')] < ^ [kolBy(Rd)] foranyt>0. (4.2) 

Proof. For fixed e > 0, we mollify uq by Uq € C°° so that E \\\uo\\Hs(j^d)\ is finite 
for any s > 0, and 

^ [|I"oIIlp(R<*) + WolBViM")] < E [llwollLp(Rd) + |wo|i3y(R^)] < oo, 

for any p = 1, 2, • • • , and Uq{:s.) — > mo(x) for a.e. x, almost surely as e 0. 

Now the same arguments as in Section 4 of Feng-Nualart [9] yield that there 
exists an ^"(-adapted stochastic process u'^ = u^{t) G C([0, oo); L^(M'')) satisfying 
almost surely that 

(i) E [\\u'{t, •)ll?/=(Rd)] < for all t > 0; 

(ii) d^^^^u'it, •) e C(M'^) for all i,j^l,...,d; 

(iii) For any (p e C^iR"^), (p > 0, and < s < t, 
{r,iu^{t,-)),ip)~{^iu'{s,-)),ip) 

\q{u^T,-)),\/p) dT + ^j\r/'{u'{T,-)){a{u^T,-)f,p) dr 
+ (77'(«(T,.))a(u(T, •)),¥>) dW{T) 

+ ej'[{7j{u^{r, .)), A^) - {v"{u^r, •))|Vzi^(r, ■)\\p)) dr 
< J\ci{uHt,-)),W^) dT+^l\r,"{u%T,-)){a{u^r,-))\^) dr 
+ f\7f{u'ir,-))a{u{r,-)),p) dW{r)+0{e), 

J s 

where the first equality in (iii) follows from the Ito formula. 

Combining the results established in Sections [2] and [3l we conclude that there 
exist a subsequence (still denoted) {M^(i,x)}£>o and a limit u(i,x) such that as 
0, 

u'^(i,x) u{t,x.) for a.e. (t,x), almost surely, 

and the limit u(t,x) satisfies (|4.2I) . Arguing as in Feng-Nualart [9], we can pass to 
the limit in the entropy inequality (iii) to conclude that the limit function u{t, x) 
is a stochastic entropy solution (cf. Definition 14. ip . Moreover, we can prove that u 
is a strong stochastic entropy solution, as defined in Definition 14.21 □ 
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Combining Theorem 4.1 with the L^-stabihty result established in Feng-Nualart 
[9], we conclude 

Theorem 4.2 (Well-posedness in L^). Suppose (|1.4p and (|1.5p hold, and that uq 
satisfies 

^ [ll"o|lip(Rd)] < oo, p = 1,2,-- - . 

(i) Existence: There exists a strong stochastic entropy solution of the balance 
law ()l.ip with initial data (II. 2p . satisfying for any t > 0, 

p=l,2,.... (4.3) 

(ii) Stability: Let u(i,x) 6e a strong stochastic entropy solution of (jl.ip wif/i 
initial data uo(x), and fei w(t, x) be a stochastic entropy solution with initial 
data Wo(x). Then, for any t > 0, 



E 



/ \u{t,x) — v{t,x.)\ dx <E / |uo(x) — i'o(x)| dx 



(4.4) 



Proof. For the n^]^LP(R'')-valued random variable uq, we can approximate uq 
by uf^{x) in as (5 ^ 0, with £;[||u;^||p + \uf^\Bv] < oo for fixed 5 > 0. Then 
Theorem 14.11 indicates that there exists a corresponding family of global strong 
entropy solutions u^{t,x) for S > 0. 

Then the L-'^-stability (contraction) result established in Feng-Nualart [9] implies 
that u^{t,x) is a Cauchy sequence in L^, which yields the strong convergence of 
M*(t,x) to u{t,x) a.e., almost surely. Since 

E [\\u'{t, OIILcM^)] < E < C, P = 1, 2, • • • , 

where C is independent of S, one can check that u{t, x) is a strong stochastic entropy 
solution, and (|43|) holds. For the stability result (|44l) . see [9]. □ 

5. Continuous dependence estimates 

The aim of this section is to establish an explicit "continuous dependence on the 
nonlinearities" estimate in the BV class. Let u{t) — u{t, x; ut) be a strong stochastic 
entropy solution of 

dtu + W ■ i{u) = a{u)dtW, u|t=o = uo- (5.1) 

Let v{t) = v{t,x;uj) be a strong stochastic entropy solution of 

dtV + V ■i{v) = aiv)dtW, v\t^o = vo. (5.2) 

We are interested in estimating E [\\u{t) — f (i)||^i] in terms of uq — vq, { ~ f, and 
a — a. Relevant continuous dependence results for deterministic conservation laws 
have been obtained in |JJfl[2], and in [3] for strongly degenerate parabolic equations; 
see also [3l [11] . 

We start with the following important lemma. 

Lemma 5.1. Suppose that (|1.3p - (|1.5p hold for the two data sets (uo,f, cr) and 
(fQjf, (t). For any fixed e > 0, let u{t) = u(t,x;u)) be the solution to the stochastic 
parabolic problem 

du+\y^-i{u)-e/^^u]dt^a{u)dW{t), u\t=o = uq. (5.3) 
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For any fixed e > 0, let v{t) = vlt,y]uj) he the solution to the stochastic parabolic 
problem 

dv -iiv) - e/:^yv\dt = a{v)dW{t), v\t=o = vo. (5.4) 
Take Q < (ps ^ 05 (x, y) e C^{W^ x W^) to be of the form: 

M^,y) - wA^W^) =■■ Jsi^W^), (5.5) 

where J(-) is a regularization kernel as in (|331) and < V G C;?°(R^). Moreover, 
given any entropy function r]{-) with 77(0) = and r]'{-) odd, introduce the associated 
entropy fluxes for u,v G K.- 



J y J y 

Then, for any t > 0, 

J J viuit,x) ~ v{t,x))<f)s{yi,y)dxdy - JJ ?7(wo(x) - wo(y))'^<5(x, y) dxdy 

^ Ills 

where 

= j I j\\u{s,x),v{s,y)) ■^i,{^)Js{^)dsdxdy, 
l''\<Ps) = jjj^ (q'^(w(s,y),w(s,x)) - q*^(u(s,x),w(s,y))) • Vy(/>5(x, y) dxdy, 
I'^'m ^{V~e-Vifll f\{u{s,x) ~vis,y)) Ay Js{^)i^{^)dsdxdy 



+ ^{^e + Viyff l\{u{s,x) - v{s,y))Js{^)Atpm^)dsdxdy 



+ {i-e) 77(u(s,x)-i;(s,y))VyJ5(x-y)-VV'(^)dsdxdy, 











X ((t(m(s,x)) - (T(u(s,y)))^05(x,y)dsdxdy. 



Proof. Subtracting (|5.4I) from (|5.3p and subsequently applying Ito's formula to 
rj(u(t) — w(t)), we obtain 



d7]{u -v)= - ri'{u - v)(Vx • i{u) - Vy • i{v)) + rj' [u - v){eA^u - eAyv) 

+ ^T]"{u-v){a{u)-a{v)y^dt (5-6) 
+ ri'{u - v){a{u) - a{v)) dW(t) 



Observe that 

V'{u - v)V^ ■ f (u) = Vx • q*'(u, v), rj'iu - v)Vy ■ {{v) = Vy • q^{v, u), 
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and thus 

-?7'(u-«)(Vx-f(w)-Vy -fH) 

= -(Vx + Vy) • q'{u, v) + Vy {q'{u, v) - q'{v, u)). 

Next, 

ri'{u — v)(^eAxU — eAyw) 

= (eAx + £Ay)r]iu - v) - r]"{u - v){s\V^u\'^ + £\Vyvf) 
= (eAx + 2^/eViVx • Vy + iAy)r]{u - v) - t]"{u - v)\y/eS/^u - Vi\/yv\^. 
Inserting the last two relations into (|5.6p . we arrive at 

d-niu -v)^ - (Vx + Vy) • q^(u, w) + Vy • (q(M, v) - q^{v, u)) 
+ (eAx + 2Ve\/iVx • Vy + eAy)?7( u — vj 

-v"iu-v)\VIV^u-ViVyv\'^ (5.7) 



1 r 

+ v){a{u) -(T{v)y 



dt 

+ ri'{u~v){a{u)-a{v))dW{t). 
We integrate (|5.7p against the test function (f)s defined in (I5.5p . yielding 



?7(u(t,x) - i;(i,x))0i(x,y)dxdy ~ j j ?/(uo(x) - -i;o(y))05(x, y) dxdy 
<ll+ll + Id + I^-\4>s) 

7j'{u{s,x) - u(s,y))(cr(-u(s,x)) - cr(w(s, y))) ^^(x, y) dW{s)dxdy, 



where 



^i-=-JJ (Vx + Vy) •q''(u,«)(/.5(x,y)dsdxdy, 

Ic JJJ^ ' - q^(w(s,y),'«(s,x)))(^i-(x,y)dsdxdy, 

Id-= [[[ (eAx + 2^/i\/iVx• Vy +£Ay)77(u(s,x)-w(s,y))(^5(x,y)dsdxdy. 



Integrating by parts gives = I^'^{(f)s), and also — I^{(f)s), since 

(Vx + Vy)0 = J.(^)(Vx + Vy)V'(^) = (^)V^(^). 

We now investigate the term Id- A calculation shows that 

(eAx + 2V^\/iVx • Vy + eAy)0(x, y) 

= (eAx + 2ViViVx • Vy + eAy)Js{^M^) 

+ Jsi^ - y)(eAx + 2Vi%/fVx • Vy + eAy)V(^) + R, 
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and 



R = 2eV^Js{^) ■ V,,V(^) + 2eVy J5(x - y) • VyV(^) 

+ 2^/iViVxJ^(^) • Vy7/;(^) + 2^eVi^yJs{^) ■ VxV^(^) 

= (2£Vx J5(^) + 2^/iViVxJ^(^) + 2ViVfVy J5(^) 

+ 2fVyJ,(^)).Vy^(^) 

= 2VyJ,(^) • VyV(^)(£ -e)= VyJsi^) ■ VV'(^)(e - e). 

Moreover, 

(£Ax + 2v^ViVx-Vy + £Ay)J5(x-y) = (V^ - \/i)^Ay J5(x - y), 
(£Ax + 2ViVfVx-Vy + £Ay)i^(^) = i(Vi+Vi)'AV'(^). 

Consequently, after integrating by parts. Id becomes P''^{(j3s)- 



□ 



Theorem 5.1 (Continuous dependence estimates). Suppose that (|1.3p ~ (|1.5p hold 
for the two data sets (uo,f, cr) and (vo,f,a). Let u(t) and v{t) be the strong sto- 
chastic entropy solutions of (|5.ip - (|5.2p . respectively, for which 



In addition, we assume that either 

u,v e L°°((0,T) X M'' X Q) for any T > 0, 



f", f -{',a- a e L" 



Then 



(i) there is a constant Ct > such that, for any < t < T with T finite. 



E / |M(i,x) - w(i,x)|V'(x)dx 



<Ct\ E 



|uo(x) - wo(x)|'0(x) dy 



Vtm 



tE [\vo\ 



BV( 



f'-f'l 



Loc + cr - cr Lo 



where the constant Ct > is independent of \uo\gy(g^d-^ and \vo\bv(R'')' 
and may grow exponentially in T. Moreover, ip = ?A(x) > is any function 
satisfying jV'l < Cq, < CqV'; which includes = e~^°^^^ and, more 

generally, ^/'(x) — 1 when |x| < R and V'(x) = e"'^"^''''"''^-' when |x| > R. 
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In particular, for any R> 0, this choice implies 



x\<R 

< Ct,r E 



\u{t, x) — v(t, x)| dx 

|mo(x) - uo(x)| dx 



Vi\\a - (7||l= 



+ tE [\vo\sy^^.)] (Ilf - f'lU^ + Ik - a|lL~) 



(ii) There is a constant Ct such that, for any < t < T < oo, 



E / |u(t,x) - ■i;(t,x)|7/;(x)dx 



<Ct\E 



|uo(x) - fo(x)|V'(x)dx 



+ V^IIV'llLi(R<i)A(CT,CT) 



+ t E [bolBy(K^)] (llf - f'lU- + a)) j , 



where V'(x) is as before and 



A(cr, a) := sup r-^ . 

Remark 5.1. If, in addition to the assumptions listed in Theorem 15. 1[ mo(x) and 
?;o(x) are periodic in x with the same period, we can "remove" ip from the above 
estimates, since integrations are then over a bounded domain. 

Proof. As the vanishing viscosity method converges (cf. Theorem 14. ip . it suffices to 
prove the result for (|5.3l) - (|5.4p with e — e. 

For /9 > 0, let 7?p : M -J> M be the function defined by (EH) -([23]). Then the 
function 



satisfies 



9„(q^(u,«)-q^(^;,u))|<^||f"IL..p, 



(5.8) 



where M2 = sup|„i<i W{u)\. 
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In view of Lemma l5. II with i = e, 



E 



< E 



rip{u{t,x.) - v{t,y))(l)s{x,y)dxdy 



E 



+ E 



E 



eE 



77p(uo(x) - ■(;o(x))05(x,y)dxdy 



c^l{u{s,^),v{s,y))-^^{^)Js{^)dsd^dy 



(^p{v{s,y),u{s,yi)) 



q|^(u(s,x),w(s,y))) • Vyc/)^ dsdxdy 



-77^'(w(s,x) - w(s,y)) 



2 

X ((7(^(5, x)) - (T(u(s,y))) (^5(x,y)dsdxdy 



77p(w(s, x) - z;(s, y)) J5(^)Ax^/'(^) ds dxdy 



Observe that 



Vy ■ (qp(w(s, y), x)) - qp(u(s, x), w(s, y))) 

= Vy« • 9„(q^(w,i') - qp(«,«))|(„,„)=(„(,,,),„(,,y)), 



and, thanks to 



= c't,(qp(w,u) ~ qp(w,u)) + 9^ (qp(u, u) - qp(v, w)) 
<|f'(«)-f'(i;)| + ^||f"|U^p. 



Hence, after an integration by parts, 



E 



(qp(«(s,y),'«(s,x)) - q^(u(s,x),w(s,y))) • Vy(t)sdsdxdy 



< 



tE [K'olw(E<<)] II^IIl~(m^) (llf - f'lU- + ^l|f"lk»p) 
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Consequently, again thanks to (|5.8p and also (|2.4p . we can write (|5.9p as 
E I I \u{t,x) - v{t,y)\(t)s{x,y)dxdy 



E 



< E 



+ E 



|uo(x) - wo(x)|05(x,y)dxdy 



1 



-77 (u(s,x) - v{s,y)) 



(5.10) 



X (cr(u(s,x)) - o-(w(s,y)))^(/)5(x,y) dsdxdy 

+ t MBViR") IIV'llL~(K<i) 



Li(K<i)P) +C(£)- 



Sending 6 and using jVV'Cx)! < Co^/'(x), we obtain 



lim 



E 



q^(u(s, x), ^(s, y)) • Wip{^)Js (x - y) ds dx dy 



<C2||f'| 



|m(s, x) — v{s, x)| ^(x)dx 



ds: 



hence, sending i5 — >■ in (|5.10p returns 



E 



\u{t, x) — v{t, x)| dx 



|uo(x) - wo(x)| V'(x)dx 



<C2||f'l 

+ E 
+ 



E 



\u{s, x) — w(s, x)| '0(x)dx 



ds 



-ri'p{u{s,x) — w(s,x))(cr(ii(s,x)) — (t(w(s, x))) ^7/;(x) ds ds 



f'-f 



Next, with our choice of r/p, it follows that 



0{p)) 



E 



< E 



-rip{u{s,x) - w(t, x))(cr(u(s,x)) - (7(i;(s, x)))^?/;(x) dsdx 



— ll«(s,x)-«(s,x)|<p(o'(u(s,x)) - (t(m(s,x)))^'0(x) ds dx 
P 



(5.11) 



=: A + B 



— l|«(s,x)-t,(s,x)|<p(<5'(u(s,x)) - o-(w(s,x)))^7/)(x) ds dx 

P 
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Clearly, 



\A\ < C3E 



* |ct('u(s,x)) - (t(m(s,x))|^ 



/o P 

< C'3||V'||Li(R'i)^lk - O'lli^/P 



■(/'(x) (is dx 



and, in view of (jl.Sp . 

\B\<Cij^E y"|u(s,x) - w(s,x)|^/;(x)dx 



ds. 



In summary, we have arrived at 



E 



\u{t, x) — v{t, x)| V'(x) dx 



|uo(x) - wo(x)| ?A(x) dx 













[/ 




Jo 





|m(s, x) — w(s, x) I ■(/;(x)(ix 



+ |1VL»(K.)£; [kolw(K^)] t (iif - f'lu^ + p) 



which implies via the Gronwall inequality that, for any t > 0, 



E 



i{t, x) — v{t, x) I ')p{'x.) dy 



< gC||f'||i,=oi^ 



|uo(x) - wo(x)| V'W d> 



(5.12) 



+ Ce^llf'lU-M [kolBy(R.)] ^ (l|f' - f'lk~ + p) 



+ IIV'llLi(R'')i|k-0-|lioo//3+ ||Vj||ii(Kd)P + f 



Choosing p = \/t||<T — (7||ioo and sending e — >■ supplies part (i). 

About part (ii), the only difference in the proof comes from the estimate of the 
A-term in (jS.llI) . which is replaced by 



\A\ < C3E 
= C3E 



* |ct(w(s,x)) - (t(u(s,x))P 



* (A(a,a)) 
P 



p\u{s, x)p 

2 



|m(s, x)|^'0(x) ds dx 



\u{s, x)p^(x) ds dx 



< C'3llV'llL==(R<i)-E [||'«||l=o(o,T;L2(R'*))] 
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With this estimate at our disposal, (|5.12|) is replaced by 



E y |-u(t,x) - w(t,x)| V(x)dx 



uo(x) - t;o(x)| -0(x) dx 



l'/'llL~(R<i)- 



Part (ii) follows by choosing p = ^/tA{a, a) and sending e ^ 0. 



□ 



Theorem 5.2 (Error estimate). Suppose (jl.3p - (|1.5p hold. Let u{t) he the strong 
stochastic entropy solutions of (15.11) . for which 



(5.13) 



and let he the solution to the parabolic problem 



du" + [Vx • f (u^) - eAxm"] = cr(ii^) dVF(i), u^|t=o = "o- 



In addition, we assume that 



either u,v & L°°((0, T) x R'' x n) for any T > 0, or f" e 



T/ien i/iere exists a constant Ct > such that, for any < t < T with T finite, 



E 



/ |u(t,x) - M''(i,x)| dx < Ct-E [|uo lay /Rd J iVe- 
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Proof. Wc proceed as in the proof of Theorem 15.11 starting off from Lemma 15.11 
with a = a,{ = {,s^e,u'^ = u, u^ = v, leading to 



E 
< E 



'(t,x) - M^(t,y) 05(x,y)dxdy 



Mis,^),u'{s,y))-Vi;i^)Jsi^)dsd^dy 



+ E 



X ((t(w''(s,x)) - cr(u''(s,y)))^(/)5(x,y)dsdxdy 

t IMsViR") IIV'llL~(Rd) 0{p) + O (||l/^||ii(Rd) p) 



{V^^Vi^E 



+ ^{V-e + VifE 



+ {e-e)E 



?7p(m^(s,x) - u^(s,y)) 

X AyJ5(^)^/;(^)dsdxdy 
r]p{u''{s,x.) - u^(s,y)) 
X Js{^)Aij{^)dsd^dy 
77p(m''(s,x) - u''{s,y)) 

X VyJ5(^) • VV^(^)dsdxdy 



= : /i +/2 +/3 + /4 + /5 +/6- 

As before, 



x) - u'{s, y)| J5(a; - 2/)V(^) dxdy 



(5.14) 



Noting that the right-hand side is independent of p, we can first send p -> in 
(15.141) . and then let ip tend to l^d, keeping in mind the LP-estimates (|4.2p . with 
the outcome that Ii, I^, I5, Iq — > 0. The resulting estimate reads 



E 



u'{t,x)-u'{t,y) Js{^)dicdy 



<Ci E 



\u{s, x) - v{s, y)| Jsi^) dx 



(5.15) 



ds + 1, 



where 



i = {^-ViyE 



\u%s,x) - {s,y)\AyJs{^) ds dyidy 
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An integration by parts, followed by application of the spatial 5y-estiinate 
(f5J3l) . yields 



\I\<C2tE[\ 



lBV(R'')\ g 

In view of this, it follows from (I5.15|) in a completely standard way that 



E 



<Ci E 
Jo 



|m^(s,x) — v^{s,x)\ dx 



ds 



+ C3E [|uolBy(R^)] + t 
Choosing 6 = \fe — \fk gives 



E 



|M^(t,x) - u^{t,yi)\dii. 



< 



CtE[ 



Vi) 



Sending e —?' concludes the proof of the theorem. 



□ 



Remark 5.2. Theorem 5.2 indicates that {u'^(i,x)} is the Cauchy sequence in 
C(0,T;L^), which directly implies its strong convergence. 

6. More General Equations 
We now discuss briefly diverse generalizations. 

First of all, as in [9], the stochastic term in (jl.ip can be replaced by the more 
general term 

/ a{u{t,x); z)dtW{t,dz), 
Jzez 

where Z is, a. metric space, cr : K. x Z — > K, W{t, dz) is a space-time Gaussian white 
noise martingale random measure with respect to a filtration {Ft] (see e.g., Walsh 
[24], Kurtz-Protter [M]) with 

E\W{t, A) n W{t, B)] = pl{A n B)t 

for measurable A,BcZ, where fi is a. (deterministic) (j-finite Borel measure on the 
metric space Z. In particular, when Z = {1, 2, . . . , rn} and fi is a counting measure 
on Z, then the stochastic term reduces to 



fc=i 



akiu{t,^))dtWk{t). 



For the spatial BV and temporal L^-continuity estimates and stability results, 
we can allow for more general fiux functions f(t,x, u) with spatial dependence, by 
combining the present methods with those in [3J [TT] . 

Next, let us discuss the case where the noise coefficient cr(x, u) has a spatial 
dependence, focusing on the stochastic balance law 



dtu + V ■ iiu) ^ a{^,u) dtW{t), 



(6.1) 
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where the noise coefficient is assumed to satisfy cr(x, 0) = and 



|(t(x, u) — cr(x, v)\ < C\u — v\ , 



yu,v e M, Vx e 



|(7(x, u) - a(y, u)\<C\x~y\\u\, Vm e M, Vx, y e 



(6.2) 



where C is a deterministic constant. 

In the previous sections, we have estabhshed the existence of a strong stochastic 
entropy solution in the multidimensional context. The proof was based on deriving 
Sy-estimates. However, as mentioned before, the BF-estimates are no longer 
available when the noise term a depends on the spatial location x. However, it 
possible to derive fractional BV estimates. For fixed £ > 0, let w'^(i,x) be the 
solution to the stochastic parabolic problem 



[Vx • f (w^) - eAxm"] dt = ct(x, u^) dW{t), u'\t=o = uo. 



(6.3) 



where we tactically assume that f , a, uq are sufficiently smooth to ensure the ex- 
istence of a regular solution [9]. Utilizing the continuous dependence framework 
(Lemma [531 which also holds when the noise term a depends on x, we will prove 
that, for any (5 > 0, 



E 



/ / \u'^{t,x + z) - u'^{t,yL- z)\ Js{z)tp{x) dyidz 

uo(x + z) — uo(x — z)| Js{z)ip{:x.) dxdz 



<CtE 
+ CtS 



(6.4) 



1 



l^ll 



Ll(R'i)y 



< t < T, 



for some finite constant C't independent of e, where Jg is a symmetric moUifier and 
■0 > is a compactly supported smooth function. In what follows, we assume that 
the cut-ofT function ip > satisfies 

|V7A(x)| <CoV'(x), |AV'(x)| < Co^(x), V'=lonKfl:={|x| <i?}. 



for some constants Co > and R > 0. One example of such a function, at least 
after an easy approximation argument, is the compactly supported function e 
p[/2,oo(j^d) defined by 



1 



V'(x) 



when |x| < R, 

when i? < |x| < i? + TT, 

when Ixl > i? + tt. 



Estimate (|6.4p can be turned into a fractional BV estimate thanks to the follow- 
ing deterministic lemma, which is related to known links between Sobolev, Besov, 
and Nikolskii fractional spaces (cf., e.g., [E]); a proof can be found in the appendix. 

Lemma 6.1. Let h : M"^ — ^ R 6e a given integrable function, r, s G (0,1), V' S 
C^(]R'^), and {Js}s^q a sequence of symmetric mollifiers, i.e., Js{x) = jdJ {^^), 
< J e C^iR), supp{J) C [-1,1], J(— ) = J(-); and S J ^ 1. 
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Suppose r < s. Then there exists a finite constant Ci — Ci{J,d,r, s) such that, 
for any S > 



/ / |/i(x + z) — /i(x — z)| J5(z)'0(x) dxdz 
<CiS'' sup |zP" / |/i(x + z)-/i(x-z)|^/'(x)dx. 

z|<(5 Jr^ 



(6.5) 



Suppose r < s. Then there exists a finite constant C2 = C2{J,d,r, s) such that 
for any (5 > 

sup / |/i(x + z) — /i(x)| ?/;(x) dx 

\z\<S JRd 

<C2(5'^ sup 5-' ( [ |/i(x + z)-/i(x-z)| J5(z)V(x)dxdz ^^-^^ 

0<<5<1 jRd jRd 

Suppose uo is say a deterministic function belonging to BV{M.'^), or more gener- 
ally to the Bcsov space B{ ^{M.'^) for 1/ E 1). 
Starting off from with 6 > 0, 



/ / X + z) — X — Z)| J5(z)?/'(x) dxdz 

Jr'' JR-i 

2 / / |uo(x + z) - Mo(x - z)| J5(z) V'(x) dxdz 

Jr<' JRd 

+ Ct (i + IIV-IIlmr-)) 
<2CtCi IIV^II^oo^Rd) sup \zf / |mo(x + z) - uo(x)| dx 

■z|<i5 JR-i 



(6.7) 



+ Ct (i + II^IIli(r.)) 

< C{T, R), 

where (|6.5p with = 5 and s > i was used to arrive at the second inequality. 
In view of (j6.6p with s = ^ and r < i, 



sup i? 
I^l<l 



/ |M^(t,x -I- z) - ^^(tjx)! V-Cx) dx 
<C2<5'' sup S-^ [ [ |u"(i,x + z) -w'^(t,x-z)| J5-(z)VXx)dxdz ^^'^^ 

0<<5<1 jRd jRd 



+ 026'' ||u'(t,-)llLi(M.)- 

Combining (16.71) with (16. 8p yields 



Theorem 6.1 (Fractional BV estimate). For fixed e > 0, let u'^ solve the stochastic 
parabolic problem ()6.3p with deterministic initial data uq belonging to the Besov 
space Bi ^{R'^) for some v G (^, 1). In addition, we assume that 

either e i°°((0, T) x M"* x n) for any T > 0, or f" e L°°. 
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Fix T > and R > 0. There exists a constant Ct,r independent of e such that, for 
any < t < T , 



sup E 

|z|<(5 



/ |u^(t,x + z) -u^(t,x)| dx 



for some r G (0, i). 



Proof of (16. 4p . We start off from Lemma [5.11 with { — i, i — e, a — a, vq — uq, 
V — u (which also holds when a depends on the spatial location): 



E 



< E 



rjp{u'{t,^)-u'{t,y))Js{^M^)d^dy 



E 



ql{u'{s,^),u'{s,y))-V^{^)Js{^)dsd^dy 



E 



E 



eE 



qpK(5, y),u%s, x)) - <^p{u%s, x),u%s, y))) • Vy(t)s ds dxdy 



-<(«^(s,x)-«-(s,y)) 



X (fT(x,M''(s,x)) - cr(y,u^(s,y)))V5(x,y)dsdxdy 



Vp{u'{s, x) - u'{s, y))Js{^)A^yj{^) ds dxdy 



(6.9) 



Finally, denoting the left-hand side of (|6.9I) by LHS and utilizing (12.41) . we have 



LHS = E 



^(i,x) - u'it,y)\ J,(^)V;(^)dxdy 



E 



\uoix) - uoiy)\ Jsi^M^) dxdy 
Since |VV'(x)| < CoV'(x), 

\h\<C \u'{s,^)~u'{s,y)\Js{^m^)d^dy 
Note that, thanks to (|5.8p and the boundedness of f", 

nu 

q^(w,u) =q^(u,«)+ / d^{cil{i,v)~-cil{v,£)) di = cil{u,v) + \u~v\0{p), 

J D 



ds. 
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SO that 

I/2I < CpE JJJ^ \u^s,^)-u^s,y)\ \VyJs{^)\^m^)dsd^dy 



CpE 



|u^(s,x) - u'{s,y)\ Js{^) |VV(^)| dsdxdy 



< Ct 



because of the estimate 



E\}\u^{t)\\ 



sup iJ/ V^yil^WKd) 
0<t<T ' 



< 00, for any T > 0, 



and we have again exploited jV-fACx)! < Cotp{x) 
Regarding I3, 



\h\<E 



P 



■l|n»(s,x)-uHs,x)|<p(o-(x, -u^(s, x)) - a{y, u^(s, x)))' 



+ E 



M. 
P 



X Js{^)^A^)dsd^dy 
^ l|u-(s,x)-«=(s,y)|<p(o-(y, ""(s, x)) - cr(y, M^(s, y)))' 



where, cf. second part of (|6.2 



* |cr(x,'u'^(s,x)) - cr(y,u^(s,x))|^ 



J5[^)^{^)dsd^dy 



< CE 



p 



u"(s,x)|V5(^)V'(^)dsdxdy 



where we have put to use the estimate 



sup E \\\u^it)\\ 

0<t<T 

Moreover, cf. first part of 



< 00 for any T > 0. 



\B\<cJ^e[JJ \u^{s,^)-u^s,y)\Js{^M^)d^dy 
Regarding I4, using |AV'(x)| < CoV'(x), we have 



I/4I <C f E 
Jo 



^{s,^)-u%s,y)\Js{^)^m^)d^dy 



ds. 



ds. 
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Summarizing, we have arrived at 



E 
< E 



«o(x) - uo{y)\ J,(^)VX^) dxdy 



C / E 

Jo 

+ Ct h 



P 

Optimizing with respect to p (take p = 0{S^^^)) and applying GronwaU's lemma 
gives 



^2 



E 



|u-(t,x) - u^{t,y)\ J,(2£-y)^(2£±y)dxdy 



< CtE 



|uo(x) -uo(y)| Jsi^)ip{^)dyidy 

+ C7T(i + ||^|lii(R.)) A o<t<r, 



(6.10) 



for some constant Ct independent of e. 

Introducing new variables, x — and z = in (|6.10p . so x = x + z and 
y = X + z, we finally obtain (dropping the tildes) (|6.4p . □ 

Combining Theorem 6.2 with the argument in Section 3, we conclude 

Theorem 6.2 (Existence and regularity). Let (|6.2p and ||f"||ioo < oo hold. 

(i) Let the initial data uq belong to the Besov space ^{W^) for some v g 
(i, 1) and 



E hllL( 



< oo, p = 1, 2, 



3.11) 



Then there exists a strong stochastic entropy solution of the balance law 
(16.11) with initial data uq such that, for fixed T > and R > 0, there exists 
a constant Ct.b. such that, for any < t < T , 



sup E 

|z|<A- 



\u(t, X + z) — u{t, x)| dx 



Kr 



< Ct,r 5^ 



for some re (0; 5) o-'^d 



3.12) 



(ii) Let Uq satisfy only (|6.1ip . Then there exists a strong stochastic entropy 
solution of the balance law (|6.ip with initial data Uq satisfying (|6.12p . 

Finally, we remark in passing that the results and techniques straightforward ex- 
tends to nonlinear stochastic balance laws with additional nonhoniogeneous terms, 
by combining with the Gronwall inequality, such as 

dtu{t, x) + V • {{u{t, x)) = a{u{t, x)) df W{t) + g{u{t, x), t, x), x e M'', t > 0, 

with initial data (|1.2p . for a large class of nonhomogeneous terms g{u,t,x.). 
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Appendix A. Proof of Lemma 16.11 
Suppose r < s, and let us prove (|6.5p as follows: 

^'"^ / |/i(x + z) - /i(x - z)| J5(z)V'(x) dzdx 
M2i±^l^^fc^ J(M)^(x) dz dx 



Sd+r 



<\\J\\ 



z\<S 



|/i(x + z) — h{x — z) 

I \d+r 



■■0(x) dz dx. 



< \\J\\L-°(m sup Z " +Z) - h{- - Z)) V-II^WRd) / 



|z|<5 

< Cj^d,r,s sup Z"'* +z) - h{- - z)) V'llil(Rd) , 

z| <(5 



|z|<5 |Z| 



where we have used the integrability of 1/ l^l'^^'^ " (since d + r — s < d). 



We continue with the proof of (|6.6p . To this end, let us introduce the modulus 
of continuity 

:= sup / |ft,(x + z) - ft,(x)| ■i/;(x) dx, 5 > 0. 

|z|<<5 JR'i 

Clearly, a;(-) is a non-decreasing function and thus 



K ^ ^uj{K)dK> j K ^ ^uj{k) dn > uj{S) j n ^ ^ dn ~ —5 ^uj{5)] 



therefore 



Set 



and note that 



uj{5) <r6'' k-'-^oj{k) dn. 
Jo 

hs{yi) := J6{y)h{x + y)dy, 



(A.l) 



/ |/l(x + z) -/l(x)|7/;(x)dx 

< \hs{x + z) - hs{x)\tl;{x.)dx+ \hs{x + z) ~ h{x + z)\ip{x.) dx (A.2) 



+ / |/i5(x)-/i(x)|7/;(x)dx, 

JR'' 

We estimate the first two terms on the right-hand side as follows: 
\hs{x) — h{x)\ ipix) dx 



2d^-d / J(2M)(/,(x + y)-/j(x))dy 



ip{x) dx 



< \\J\\ 



/ / |/i(x + y)-/i(x)|V(x)dxdy 

J|y|<| Jm" 
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and, similarly, 

/ |/i5(x + z)-/),(x + z)|V'(x)dx 

,5-'*/ / |/i(x + z + y)-/i(x + z)|V(x)dxrfy 



lz,~(M) ^ 



\y\<i 

' [ f |/i(x + y)-/i(x)|^(x-z)dxdy 

<C6-'^ f f |/i(x + y)-/i(x)|V(x)dxdy + /i((5), 

^|y|<| Jr'^ 



where, for (5 > 0, 



7i(^) := S-" sup / / |/i(x + y) - /i(x)| |^(x) - V(x - z)| rfxdy 

|z|<| "'|y|<'5 "'R'* 

< (5C||VV'|lio.(Rd) ||/i|lLi(Md) lo<5<l(<5) 
+ IIV'llLoo(Rd) ||/l|Li(Rd) l5>l(5). 

For each z e R"^ and x e M"^, 

hs{x + z) - hsix) = I Vhs{x + £z) ■ zd£. 
Jo 

Observe that for each a; e R"^, 

V/i5(x) = / V J. (y) (/i(x + y) - /i(x)) rfy. 

by the symmetry of the moUificr. Thus, with |z| < 6, 
/ \hs{x + z) - hs{x)\^li{x)dx 

= [ [ V/i5(x + £z) • zd£ V(x)rfx 
Jr<i Jo 

<C6-'^ sup / / |/i(x + £z + y) -/i(x + ^z)|^(x)dxdy 

|z|<5,£e[0,l] J|y|<| JR-i 

= CS-'^ sup / / |/i(x + y)-/i(x)|V(x-fe)dxdy 

|z|<5,^e[o,i] J\v\<^ Js.'i 

KCS-" f f |/i(x + y)-/i(x)|V(x)dxdy + /2(5), 

J\y\<^ Jr.'' 

where /2(^) denotes the expression 

C6-'' sup / / |/i(x + y)-/i(x)||^(x)-V(x-^z)| rfxdy, 

\z\<S,i£lO,l] J\y\<^ Jm-i 

and 

HS) < '^C'||VV'|lz,oo(Rd) ||/t|lLi(Rd) lo<5<i('5) 
+ C ll^llLoo(Rd) ||/l||Li(Rd) l5>l(5), 

cf. the term Ii{S). 
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In view of the estimates derived above, taking the supremum in (jA.2|) over |z| < S, 
we have estabhshed 

u{5)<C5~'^ I I |/i(x + y)-;i(x)|?/'(x)dxdy 

+ C Mlur^) [s 1o<5<i{6) + l5>i{6) 

Multiplying this by 5~^~^ and integrating yields (replacing y by z) 



<C f 5-'-^-'^ f f \h{x + z)-h{x)\tlj{x)dxdzd5 (A3) 

1 poo 

r-1 



+ C|l'i|lLi(M^) (^y 6'"d5 + J 6-'-'dSj=:A + B, 
where the integrals on the last line are bounded since r G (0, 1): 

B <Cr ll^llLi(Kd) ■ 

Since x — 5^"^(t') ^"^^ remembering r < s, 
A<Cj [ 5-'-^-'^ ( [ |/i(x + z)-/i(x)| j(^) V(x)dxdzd(5 

Jo J|z|<|JR<i 

<Cj [ S-'^S''-'^-^ [ [ \h{x + z) - h{x)\Js{2)ilj{x)dxdzdS 

f' ITT^'^^] sup / / |/i(x + z)-/i(x)| J5(z)dxdz 

\Jo 0'-+'^ J o<S<l Jri jRd 



< 



Cj^r,s sup 5 I I |/i(x + z) — /i(x)| J5(z) -(/'(x) dxdz, 

0<5<1 jRd jRd 



where Cj^r,s = C'j 7^7- 

Consequently, from (|A.ip and (jA.3[) it follows that for any (5 > 0, 

sup / |/i(x + z) — /i(x)| V'(x) dx 

|z|<<5 JR'' 

<C(5''sup (J-'* / / |/i(x + z)-/i(x)| J5(z)V'(x)dxdz (^-"^^ 

0<5<1 jRd jR^i 

+ CJ'' ll^llil(Rci) , 

for some finite constant C . 
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Finally, observe that 

/ / |ft.(x + z) - /i(x)| J5(z)'0(x) dxdz 

= / / |/i(x + z) - /i(x - z)| J5(2z)V'(x - z) dxdz 

= o'd / / l^(x + z) -/i(x-z)|,/5/2(z)V'(x-z)dxdz 

/ |/1(X + Z) -;i(x-z)| J5/2(z)V'(x)dxdz+/3('^), 

where hiS) denotes the expression 

/ / |ft,(x + z) -ft,(x-z)| J5/2(z)|-0(x) -'i/;(x-z)| dxdz. 

^ Jr" JR-i 

As with h{6), hiS) < C ||/i||ii(R.) ((5 lo<5<i(<^) + l5>i(<5)) , 



and as a result 



sup S " / |/i(x + z) - /i(x)| J5(z)?/'(x) dxdz 

0<<5<1 jRd jRd 

<C sup (5"" / / |/i(x + z) - /i(x - z)| Ji-/2(z)V-'(x) dxdz 

0<5<1 jRd jRd 



■c\\hh,. 



We can therefore replace (jA.4p by 

sup / |ft,(x + z) — /i(x)| ■0(x) dx 

\z\<S ■/Rd 

<C(5''sup J-' / / |/i(x + z)-/i(x)| Ji(z)V'(x)dxdz 

0<<5<1 jRd Jk'' ^ 

for some finite constant C, which implies (j6.6l) . 
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